Sensitivity and back-action in charge qubit measurements by a strongly coupled 
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We consider charge-qubit monitoring (continuous-in-time weak measurement) by a single-electron 
transistor (SET) operating in the sequential-tunneling regime. We show that commonly used master 
equations for this regime are not of the Lindblad form that is necessary and sufficient for guaranteeing 
valid physical states. In this paper we derive a Lindblad-form master equation and a corresponding 
quantum trajectory model for continuous measurement of the charge qubit by a SET. Our approach 
requires that the SET-qubit coupling be strong compared to the SET tunnelling rates. We present 
an analysis of the quality of the qubit measurement in this model (sensitivity versus back-action). 
Typically, the strong coupling when the SET island is occupied causes back-action on the qubit 
beyond the quantum back-action necessary for its sensitivity, and hence the conditioned qubit state 
is mixed. However, in one strongly coupled, asymmetric regime, the SET can approach the limit of 
an ideal detector with an almost pure conditioned state. We also quantify the quality of the SET 
using more traditional concepts such as the measurement time and decoherence time, which we have 
generalized so as to treat the strongly responding regime. 

PACS numbers: 73.23.Hk, 03.67.Lx 
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I. INTRODUCTION 

The single-electron transistor— (SET) has been 
suggested^ as a device to measure the state of charge 
qubits, one proposal for the fundamental elements of 
a quantum computer— In solid-state systems it is well 
known that realistic measurement devices such as the 
SET and quantum point contact (QPC) do not perform 
instantaneous measurements. Rather, it is necessary to 
treat the measurement as a sequence of weak measure- 
ments, resulting in a continuous process in time. 

During such a continuous measurement, the condi- 
tioned qubit state is gradually projected into an eigen- 
state in the measurement basis. Here the conditioned 
state is one that an experimenter would calculate know- 
ing the output of the detector. On average (that is, ig- 
noring this information) the effect of the measurement 
is simply to decohere (that is, remove the coherences 
of) the qubit in the measurement basis. This quantum 
back-action is unavoidable because the coherences be- 
tween the eigenstates (which indicate that the qubit is 
in a superposition of these eigenstates) must vanish for 
one or the other eigenstate to be realized. For an ideali 
or quantum-limited^ detector, the decoherence is equal to 
the minimum back-action allowed by quantum mechanics 
given the information obtained from the detector. Conse- 
quently, a qubit state conditioned upon the output of an 
ideal detector would remain pure if it began pure. This 
is the case for the usual model of the QPC— *2*i2i 

In this paper we maintain the purity of the conditioned 
state as the ultimate quantifier of the quality of the mea- 
surement. However, it has been common to use a differ- 
ent quantifier, namely a comparison between the decoher- 
ence rate Td and the measurement rate r~ x . These are 
the characteristic rates of the decoherence and the mea- 



surement (gradual projection) processes discussed above. 
(It is assumed that there is only one characteristic rate). 
For so-called "weakly responding"— detectors, the ratio 
of these rates is limited by£> 
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with equality implying an ideal detector. While not 
the main focus of this paper, we show that in general 
the detector quality cannot be captured by such a sim- 
ple ratio, which is sometimes referred to as the detector 
efficiency^ For example, for a QPC it is easy to verify 
that TdT m varies between 1 (in the strong-response limit) 
and 1/2 (in the weak-response limit). This is discussed 
in Appendix [XI 

To date, theoretical treatments&iSii&iii^ of charge- 
qubit continuous measurements by SETs operating in the 
sequential-tunneling mode have focusscd on the weakly 
responding detector. This is where the variation in the 
detector's output that depends on the qubit state is small 
compared to the average part. We will quantify this later. 

In this paper we begin by showing that the master 
equations derived for the SET-monitored charge qubit in 
Refs. H E Q [H and 03 are not of the Lindbladii 
form that is necessary and sufficient for guaranteeing 
valid physical states. Specifically, we show that for short 
times the qubit state matrix (density matrix) may be- 
come non-positive. Positivity of the state matrix is a 
fundamental requirement for a model to represent a real 
physical system. Consequently, their evolution (the qubit 
density matrix) can violate positivity and hence cannot 
possibly represent the state of a real physical system. 
As well as our Lindblad-form master equation, we also 
present the associated equation for the conditioned qubit 
state, known as a quantum trajectory equationi^Si^ or 
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quantum filtering equationmi^ To derive our equations 
requires assuming strong SET-qubit coupling, as we will 
discuss. 

Using our model we show that the strongly coupled 
SET can, in the regime of strong response, approach 
operation at the quantum limit where the purity of the 
qubit state is preserved in a conditional measurement. In 
this regime the SET also approaches the quantum limit 
of efficiency for a charge qubit measurement given by 
Eq. p. 

The paper is organized as follows. The next section 
contains a brief discussion of the SET and charge qubit 
system. In Sec. IHII we discuss previous models for SET- 
monitorcd charge qubits and show that the master equa- 
tion in this case^ is not of the Lindblad form. In Sec. IIVI 
we present new general definitions for the decoherence 
and measurement times that are valid for arbitrary cou- 
pling strength and response. We present our quantum 
trajectory model that corresponds to charge qubit mea- 
surement by a strongly coupled SET in Sec. In the 
same section we present quantitative analyses of our mas- 
ter equations (conditional and unconditional), focussing 
on comparing sensitivity and back-action in the charge 
qubit measurement. The paper is concluded in Sec. IVII 



II. SYSTEM 

In this section wc outline a quite general model for a 
SET monitoring a charge qubit, independent of coupling 
strength at this stage. For specificity, we consider the 
double quantum dot (DQD) charge qubit^ but expect 
the results to hold for other realizations of a charge qubit. 
A schematic of the DQD-SET system is shown in Fig. ^ 
We discuss this figure below. 

For measurement purposes, sequential tunneling pro- 
cesses are the most important as they yield the largest 
SET currents^ Thus, we assume that the dominant form 
of transport through the SET is sequential tunneling, 
with higher-order processes such as cotunncling exponen- 
tially suppressed. In this regime the so-called orthodox 
theoryiS applies. This requires the SET tunnel junc- 
tion resistances to be greater than the resistance quan- 
tum h/e 2 w 26kn£22. and the SET to be biased 'near' 
a charge degeneracy points Higher-order processes be- 
come more important away from a charge degeneracy 
point, where sequential tunneling events become less en- 
ergetically favorable. There is a subtlety here for the 
strongly responding SET that we discuss when we intro- 
duce our model in Sec. 

For low temperature k^T and not-too-high bias volt- 
age eV = /il — Mr (relative to the SET island charging 
energy Ec = e 2 /2C, where C is the island's total ca- 
pacitance), the SET island has only two possible charge 
configurations, H = and H = 1. With these conditions 
satisfied, the gate voltage (which is modified by the qubit 
electron) controls the amount of current flow through the 
SET. This is what enables the SET to measure the qubit 



state. For convenience, we assume that electron trans- 
port occurs only in the source-drain direction due to the 
finite SET bias voltage eV (see Fig. QJ. 




DETECTOR 

FIG. 1: Schematic of a double quantum dot (DQD) qubit and 
coupled SET. A single excess electron is shared by the DQDs. 
When the near (far) QD is occupied, the electron tunneling 
rates through the SET junctions j = L, R are denoted T'j 
(Tj). The charging energy of the SET island is increased by 
X (h = 1) when the near QD is occupied. 

The SET source and drain leads (reservoirs) are bi- 
ased to Fermi levels /zl and /ir, respectively. Associated 
with each of the qubit states is a tunneling process onto, 
and a tunneling process off, the SET island. These pro- 
cesses occur at rates denoted Tl, T' l (source to island) 
and Tr, (island to drain). A dash indicates that the 
qubit electron is localized in the nearby QD, which we 
refer to as the "target" i21 The average current through 
the SET is I' = eT' h T' K / (T' L + T' R ) when the target dot 
is occupied, and I = eT^T^/ (Fl + Tr) when the tar- 
get dot is unoccupied. For later convenience, we now 
define the averages Tl.r = (r^ R , + Tl,r) /2 and differ- 
ences Ar L ,R = r^ R - r L ,R. 

The Hamiltonian for a charge qubit with energy asym- 
metry e and tunnel-coupling strength f2 ( sce Fig. ^) is 
Hqb = {s&z + ^offi) /2, where a x ^ z are Pauli matrices. 
We choose units such that K = 1. For e = 0, and in 
the absence of measurement, the qubit electron tunnels 
between the DQDs coherently (that is, existing in su- 
perpositions of the two localized states within the dots). 
For characterizing the measurement quality in Sec. IIVI 
it is simpler (and commonplace&iSiiiii^) to consider the 
limit of n ~~ * 0, which wc do. However, our master 
equations (both conditional and unconditional) are valid 
for non-zero Slo- When both the SET island and tar- 
get dot are occupied, the system energy increases by the 
mutual charging energy, x- This is described by the cou- 
pling Hamiltonian H x = \^ ® ^ where H = b^b and 
h = (a z + l)/2 are the respective occupation number op- 
erators of the SET island and target dot. For simplicity 
we do not show the standard Hamiltonians for the source 
and drain leads, since we trace the leads out of the total 
state matrix in order to consider only the SET and qubit 
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in the Born-Markov approximation. The Hamiltonian for 
the SET-DQD system (in the charge basis) is 

H = 1 ® H qh + H x , (2) 

where 1 is the 2x2 identity matrix. It is simple to 
diagonalize the qubit Hamiltonian to find the computa- 
tional/logical basis£ However, we are interested in the 
measurement process which, for the purposes of this pa- 
per, is performed in the charge basis. 



III. PREVIOUS WORK 

In this section we discuss previous work on the contin- 
uous monitoring of a charge qubit by a SET. We demon- 
strate that the master equation of Ref. y (and those 
of Refs. [H and HHIT^ since they are equivalent) can 
produce non-physical results and therefore is not of the 
Lindblad form. This master equation was derived inde- 
pendent of the coupling strength \, but the SET quality 
analysis always assumed weak coupling. First, we clarify 
the weak- response and weak-coupling assumptions. 

Qualitatively, a weakly responding detector is one in 
which the variation in the detector's output that depends 
on the qubit state is small compared to the average part. 
Korotko\ii2i2£ expresses this quantitatively in terms of 
the current output as 

|A/| « /, (3) 

where AI = I— I' and /= (/' + I) /2. Note that | AJ| < 
27, so that strong response is |A/| ss 21. 

The strength of the SET-qubit coupling is quantified 
by x ( see Sec. [HJl. When analyzing the quality of the 
SET, the authors of Refs. H [H Hi and El make the 
weak coupling assumption 

X « f L + f R . (4) 

In Ref. Q this assumption was made "for definiteness" . 
It leads to a major simplification in that the SET island 
charge fluctuation spectrum is 'white' over all frequen- 
cies of interest. The variations in the SET tunneling 
rates are proportional to ^^i^i^i^i^ S o that weak cou- 
pling implies weak response. However, a strongly coupled 
detector can have strong or weak response. 



A. Unconditional master equation 

The unconditional (or nonselective, or ensemble aver- 
age) master equation of Ref. is obtained by averaging 
over, or ignoring, the measurement record (the number of 
electrons that have tunneled into the SET drain) . This is 
achieved by setting k = in equations (20)-(27) of that 
paper. To perform a quantitative analysis of these equa- 
tions, it is imagined that the Josephson coupling energy 



can be "turned off" during the measurement (the equiv- 
alent Hamiltonian parameter from Sec. ^ is fio). Using 
our notations for clarity, the resulting unconditional mas- 
ter equation [equations (31) and (35) of that paper] can 
be written as 

Pn = -lLP?i+lRpii , (5a) 

f>22 = ~ r LP22 + r RP22 . ( 5b ) 
Pl2 = i£Pl2 - r L P?2 + f RPL , (5C) 

Pu = T'lPu - I>n , (6a) 

P22 = r LP22 - r RP22 , ( 6b ) 

Pu = i(e + X) P12 ~ f Rp\ 2 + ^Pi2 ■ (6c) 

Here pf^ is the element in the ith row and jth column of 
the qubit state matrix conditioned by the SET state H = 
0, 1 (the occupation number of the relevant SET island 
level). The states are normalized so that P N = Tr [p N ] is 
the probability that the SET state is H. Thus the average 
qubit state is p° + p 1 . 

We now demonstrate that Eqs. JSJ and © can produce 
a non-positive, and therefore non-physical, qubit state 
matrix. We begin by noting that the following inequality 
is satisfied by all positive 2x2 matrices. 

1 1 2 

P11P22 > |Pia| ■ ( 7 ) 

Now consider the short-time solution for p° [Eqs. JSJ] 
when the SET is initially occupied (p° = at t = 0) and 
the qubit is in the following superposition state: p\\ = 
P22 = P12 = h- After an infinitesimally short time 5t, we 
have 

P ° u « ±r' K st , 

p° 22 a ir R< 5i , (8) 
p\ 2 a ±r R Ji . 

Substituting these into Eq. J7J), canceling factors of 
(St) /4 on both sides, and rearranging gives 

|Ar R | 2 < , (9) 

which is false for all values of T R 7^ Tr . The case of 
T R = Tr corresponds to there being no response from 
the SET for a change in the qubit state, i.e. no measure- 
ment is being performed^ The result © shows that the 
inequality J7J is not satisfied, positivity is not preserved, 
and the master equation of Ref. is not of the Lindblad 
form. This is also true of the equivalent master equations 
presented in Refs. flU IT3 . ITU andll^. Numerical results 
show that the violation of Eq. (7J occurs only for short 
times (compared to the characteristic qubit decoherence 
and measurement times). However, since the qubit state 
matrix for later times is determined by the (non-positive) 
state matrix for short times, one cannot be confident in 
the ability of the state matrix at later times to describe 
a real physical system. Therefore, we maintain that the 
master equations of Refs. 0, 0, EL andlul cannot be 
relied upon at any time. 
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IV. QUANTIFYING THE SENSITIVITY AND 
BACK-ACTION 

In the previous section we showed that previous mas- 
ter equations for the SET-monitored charge qubit are not 
of the Lindblad form. In Sec. |3 we derive a Lindblad- 
form master equation (and "unravel"— it into a condi- 
tional master equation), which requires strong coupling 
between the SET and the qubit. In the strong coupling 
regime it is natural for the detector also to be strongly re- 
sponding since, as mentioned, the coupling energy largely 
determines the detector's response. The more traditional 
measures of sensitivity and back-action in this type of 
measurement aroiiSiiii^ii^ the measurement time r m 
and decoherence rate Td, respectively. Previous defini- 
tions of these quantities assumed weak detector response, 
cither explicitly or implicitly. Therefore it is necessary to 
use more general definitions, which we give in this sec- 
tion. 



A. Decoherence Time 

Previously the decoherence rate Td was defined as the 
rate at which the coherences in the qubit decay according 
to the unconditional master equation. For charge qubit 
measurements by a QPC, this description is unambigu- 
ous: there is one decay rate for the coherences which can 
be obtained^ by inspection of the master equation. This 
is not the case for the SET because there are two qubit 
state matrices, namely p° and p 1 , conditioned by the SET 
island state. Our master equation (|19fl gives a rate of de- 
cay for each of these conditional qubit coherences, neither 
of which alone accurately represents the average decoher- 
ence of the qubit state p = p°+p 1 . In Refs. 0,130! an d 
fl5L the qubit decoherence rate is defined as the slowest 
of these two rates under the premise that the qubit state 
cannot collapse any faster. This therefore represents an 
over-estimate of the quality of the measurement device, 
as it minimizes Eq. 

In an effort to more precisely quantify the decoher- 
ence, wc find a characteristic qubit decoherence time by 
borrowing a technique from the field of quantum optics. 
In a typical quantum optics scenario, where the detec- 
tor signal consists of absorption of photons emitted by a 
quantum system, the coherence of the source is defined by 
the Glauber coherence function— For a two-level quan- 
tum system (such as the two-level atom) this is given 
in terms of the steady-state, two-time correlation func- 
tion g{r) = (a-t^erO) . Here a = |0)(1| is the lowering 
operator for the atom (in the Heisenberg picture). The 
expression involves the lowering operator because the de- 
tection involves absorption. 

In the solid-state context where the detection does not 
"lower" the DQD state it is more appropriate to consider 
an expression that treats |0) and |1) symmetrically. The 



obvious expression is 

5 (r) = <at(r)a0> ss + (a(r)at0>; . (10) 

In the Schrodinger picture, this coherence function can 
be expressed as 

g(r) = Tr[*V( T )] + Tr [ap+(r)}* 

= P2i{t) + [pUt)}* , (11) 

where p(r) = exp (£ p t) p(0), p + {0) = ^ p ss , and p~(0) = 
ap ss . The Liouvillian C p is the qubit time evolution su- 
peroperator [analogous to C in Eq. I|19|)]. 

The coherence function thus defined has the property 
that it is normalized in the sense that g(0) = 1 and 
linv^oo g(r) = 0. Therefore wc can now define the de- 
coherence time (which means the same as the coherence 
time) as2& 

poo 

T d = / 9 (T)dT. (12) 

Jo 

The decoherence rate is simply Td = r^ 1 . For an expo- 
nentially decaying coherence it is easy to verify that Td 
is the exponential decay coefficient. For example, for the 
QPC master equation of Ref. Hoi this coherence function 
technique yields the same decoherence rate as the QPC 
models of Refs. 0,111 andlTol (see Appendix \K$. 

B. Measurement Time 

Previously^iSii^ii^ the measurement time r m was de- 
fined as follows. Keeping track of the number of elec- 
trons iVfi(t) that have tunneled into the drain during 
the measurement process (while ignoring the informa- 
tion contained by electrons tunneling from the source 
onto the island) allows tracking of the dynamics of the 
probability distribution for this number, P(A^r, t). Gaus- 
sian statistics were assumed for TVr, thereby implicitly 
assuming weak detector response (many SET tunnel- 
ing events occur before the qubit states are distinguish- 
able). After starting the measurement as a delta func- 
tion at Nr = 0, P(NR,t) displays two peaks (corre- 
sponding to the two qubit states) which drift linearly 
in time to positive values of N-r with velocities T = I / e 
and r' = I' /e, respectively. The widths of the peaks 
broaden due to counting statistics, growing as y/ fTt and 
y/f'T't, respectively— In Refs. and these are 
given with an extra factor of y/2 as y/2fTt and y/2fT't. 
Here the Fano factors are / = (rf, + r|) / (r L + T R ) 2 
and /' = (rj 2 + rg) / (r^ + r^) 2 . The peaks emerge 
from the broadened distribution when the separation of 
the distributions is larger than their widths, |r — r'| t > 

(vTH+zm). 

Fundamentally, the measurement time relates to dis- 
tinguishing the qubit states using information gained 
from the detector. This obviously involves the condi- 
tioned state of the qubit. Since the quantum trajectory 
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equation is by definition the optimal way to process the 
detector current to get information about the qubit, it 
should be used to define the measurement time. Our 
quantum trajectory equation is given in Sec. [V] 

In the measurement basis we can, quite generally, de- 
fine an uncertainty function for the qubit state as the 
conditional variance in a 7 : 



V(t) 



E 



= 1-E[z c 2 (0 



(13a) 
(13b) 



Here the subscript c denotes conditional variables, E[ ■ ] 
denotes a classical expectation value, and (•) denotes a 
quantum average (we use z as shorthand for (<7 Z ))- The 
initial state is the unconditional steady state, that is, the 
long-time limit of the master equation given in Sec. [V] 

Provided that the qubit Hamiltonian commutes with 
the measured observable, i.e. a z , the measurement will 
not disturb the qubit populations (but will necessarily 
affect the qubit coherence). Thus, for a non-disturbing 
measurement we require fio = 0. In this case, V(t) will 
have the properties that V(0) = 1, and linv-joo V(t) = 
0. That is, we start off with no knowledge of a z and 
end up with complete knowledge. Thus, in analogy with 
Eq. i|12|) , it can be used to give a general definition of the 
measurement time as 



V(r)dT. 



(14) 



However, because V(t) arises from stochastic evolution, 
it is not generally possible to obtain an analytical expres- 
sion for V(t), unlike g(r), which arises from the deter- 
ministic evolution of the master equation. If V(t) decays 
exponentially with time, then r m can be evaluated as 



T,„ = — 



dV(t) 



dt 



(15) 



This eq uation (expressed differently) was first used in 
Ref. |23- It has the advantage that it can be evaluated 
analytically. (In Appendix ^ we report this calculation 
for the QPC.) As we will see later, for the SET in a wide 
range of parameter regimes, the exponential approxima- 
tion for V(t) is a good one, so that Eq. (|T3|) can be used. 



V. STRONGLY COUPLED SET 

Using the quantum trajectory approach developed in 
quantum optics^ we describe both the conditional^^ 
and better-known unconditional^ time evolution of the 
charge qubit state whilst undergoing continuous mea- 
surement by a SET. In order to consider the tunneling 
processes through the SET as distinguishable, that is, in 
order to ignore the finite widths Tl,r and T' h R of these 
tunneling processes, requires^! 



This assumption of a strongly coupled SET is critical 
to our ability to derive a master equation using our ap- 
proach. As mentioned in Sec. [HI the qubit energy asym- 
metry is shifted by x when both the SET island and the 
target dot are occupied. For the strongly coupled SET, 
this causes a large shift (oc x) m the qubit Rabi frequency 
(f2 = ^/ S1q + e 2 ). One could conceivably posit that this 
situation cannot be considered as measurement as it has 
such a large effect on the qubit state. However, because 
this energy shift is in the same basis as the measurement, 
it affects only the qubit coherence and not the popula- 
tions. We analyze this below. 

Note that the time to determine the qubit state (the 
measurement time) is related to the time it takes to dis- 
tinguish between the two corresponding currents through 
the SET. Intuitively, larger separation between these 
two currents decreases the measurement time. Indeed, 
this leads to the prediction that the strongly respond- 
ing SET would enable faster charge qubit measurement 
than its weakly responding counterpart. Qualitatively, 
the largest separation of these currents occurs when tun- 
neling through the SET switches between on and off 
for the two qubit states. When the SET is operated 
near this so-called "switching point" , thermal fluctua- 
tions and higher-order processes such as cotunncling be- 
come the leading contribution to the current £ As these 
processes yield very small currents^ they may well be 
negligible in a laboratory setting. Nevertheless, we as- 
sume that the SET conducts for both qubit states to 
ensure that higher-order processes are negligible. Figure 
13 shows this assumption schematically on a hypothetical 
SET conductance-oscillation plot. The qubit state shifts 
the gate voltage by an amount that changes the SET 
conductance significantly for the case of strong response, 
while remaining always in a conducting state. 

We represent the combined state of the SET island and 
qubit by the composite state matrix 



G{t) = |0)(0| S ET®P°W + |1)(1 



SET 1 



>p\t), (17) 



x » r L + r r 



(16) 



where p(t) = Tr [G(t)\ = p°(t) + p l {t) is the qubit state 
matrix, the superscript again referring to the SET island 
occupation (N = or 1). We reiterate that, although we 
consider the state of the SET island in a composite den- 
sity matrix with the qubit, no charge superposition states 
exist on the SET island. We present master equations for 
G(t) for compactness, but could equally well present cou- 
pled maste r eq uati ons for th e qubit states p° (t) and p 1 (t) 
as in Refs. 1 1 2l and ll4ll5ll6l for example. 



A. Unconditional master equation 

Here we present our unconditional master equation for 
the state matrix G{t) of the system (qubit plus SET is- 
land). We will first outline the details of its derivation, 
which closely follows the derivation in the appendix of 
Ref. |2?]. We first consider the state matrix W(t) for the 
combined system (qubit, SET and leads). The leads are 
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V, (arb.) 



FIG. 2: Hypothetical SET conductance vs SET gate voltage. 
The two dots show the SET conductance corresponding to 
the two qubit eigenstates for the case of strong response. 



treated as perfect Fermi thermal reservoirs with very fast 
relaxation constants. Each lead then remains in ther- 
mal equilibrium at its respective chemical potential. We 
define a time interval At that is long compared to the 
relaxation time of the leads, but short (infinitesimal) 
compared to the dynamics of the SET island and qubit. 
Transforming to an interaction picture with respect to 
H [Eq. J2Jl] leaves only tunnel-coupling terms (between 
the leads and the SET island) in the interaction picture 
Hamiltonian iJj. The change in W(t) from time t to 
t + At to second order in the tunnel-coupling energy is 
then given by 



W(t + At) = W(t)-iAt H u W{t) 

rt+At 

hm, m'),w{t') 



-At 



d£ 



(18) 



We now make a Markov approximation for the instan- 
taneous relaxation of the source (L) and drain (R) leads 
so that W(t) = G(t) <g pl <g pr, and we find the result- 
ing evolution equation for G(t) by tracing over the leads. 
Assuming that all tunneling occurs from source to drain 
(fcrjT <C eV) and moving back to the Schrodinger picture, 
we obtain the unconditional master equation for G(t) as 



evolution the measurement-induced dccohercncc. It is 
defined as^S 

V [x] G = J [x] G - A [x] G , (20) 

where J [x]G = xGx\ and A [x] G=\ (x^xG + Gx^x). 
Goani& expressed the master equation of Eqs. (0 JHl in 
a similar form to Eq. I|19|) . The crucial differences, i.e. 
the non-Lindblad terms in Eqs. are shown in 

Appendix [B] 

The unconditional master equation (|19|l describes the 
system evolution when the measurement result is ignored 
or averaged over. When quantifying the SET measure- 
ment quality, we assume a non-disturbing measurement 
as noted in the paragraph preceding Eq. I|14|) . That is, 
we assume = (as was assumed in previous quan- 
titative analyses&i2ii£ii& of the SET measurement qual- 
ity). In order for our master equation to have a unique 
steady-state solution, we require nonzero f2o- Thus, for 
numerical simulations we use Qq = O.OITl. Note how- 
ever, that Eq. i|19J) is valid for any fio "C X- The numeri- 
cal solution of the master equation (|19|l gives oscillating, 
decaying off-diagonal elements of the qubit state matrix. 
We quantify this by the coherence function g(r) given in 
Eq. For SET monitoring of the qubit, C acts on the 
qubit-plus-SET state G, so p must be replaced by G in 
the coherence function of Eq. The resulting g(r) is 

shown in Fig. [21 

The four subplots in Fig. [21 correspond to increasing 
values of SET response |A/| //, which is bounded above 
by 2. The strength of the SET response was modified by 
varying T' L — the slowest rate of tunneling onto the SET 
island. The initial qubit state is completely mixed: 



Pss = -1, 



(21) 



and the SET island is initially in the corresponding 
steady-state, which is represented by the island occupa- 
tion probability P 1 = Tr [p 1 ] (with P° = 1 - P 1 ) as 



G 



H, G 



T' h V 



G + Y' K V 



b (g> h 



+T L V 
CG, 



S f (g) (i - h) 



G + T R V 



6® (i - h) 



G 
G 

(19) 



which is explicitly of the Lindblad form, and we have 
omitted time arguments of G for brevity. We remind the 
reader that b = |0)(1|set, n = (a z + 1) /2 is the tar- 
get dot occupation operator, and H is given by Eq. J5J). 
In deriving Eq. i|19f> . our assumption of strong SET- 
qubit coupling (x S> Tl + Tr) is the main differ- 
ence between our assumptions and those of the previ- 
ous modelsAi2iiiiiS4i& This assumption is essential to our 
ability to derive our master equation, as it allows each 
lead to be treated as two independent baths that differ in 
energy by h\- The Lindblad superoperator V in Eq. I|19(l 
represents the dissipative, irreversible part of the qubit 



Pi 



Ft 



(22) 



The strongest SET response [in Fig. [21 (d)] corresponds 
to 



o < r' « r L ~ r R ~ r' 



r ■ 



(23) 



where we assume that T' L is large enough in order to 
ignore higher-order processes such as cotunneling. 

The rapid oscillations in Fig. [21 are due to the strong 
SET-qubit coupling given in Eq. H16J) . Thus, these os- 
cillations are very fast compared to the tunneling rates. 
As a consequence, even if an individual tunneling event 
could be detected, it would be impossible to know when 
it occurred with sufficient accuracy to allow the phase 
of the oscillation to be known. This is a direct result 
of the strong coupling assumption (|16|) required for our 
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Lindblad-form master equation. In the observable qubit 
state, the rapid oscillations are effectively averaged over. 
This produces the observable dccoherence given by the 
bold line in Fig. [3] 

The observable state can be obtained from G(t) by 
moving to an interaction picture with respect to H x , 
unique from the interaction picture mentioned in the 
derivation of Eq. I(19|) . and making a rotating wave ap- 
proximation by dropping all rapidly oscillating terms 
(those involving x). The result is the following master 
equation 



which gives the qubit dccoherence time of 



G = -i 



T' h V P®h G 



T' K V b®h G 



g + r»v 



6® (1 -h) 



= CG, 

where the observable state G is given by 



G= |0)(0|SET®P° + |1)(1|SET®P 1 , 



G 

(24) 
(25) 

with p° = p° and p 1 = diag (/O^p^)- The Hamilto- 
nian in this interaction picture (after the rotating wave 
approximation) is 



Hi = |0)(0| S ET®flqb + |l)(l|sET' 



e „ 
2 



(26) 



Thus we see that in the interaction picture (which reflects 
the observable effects of the actual dynamics), no qubit 
coherence exists when the SET is occupied. That is, the 
off-diagonal elements of p 1 are equal to zero. 




FIG. 3: Decoherence of a charge qubit undergoing continuous 
measurement by a strongly responding SET. The rapid oscil- 
lations are too fast to be observable — the thick line shows 
the observable result. From (a) to (d): \AI\/I = 1/10, 1/3, 
1, and 3/2. Other parameters were S7o = O.OIFl, x = 50Fl, 

r R = i.ir L , r' R = i.2r L . 

Using the observable state (|25|l we obtain the following 
analytical expression for the observable coherence func- 
tion 



ffob s ("r) = P° exp (-r L r) 



(27) 



Td 



gobs(T)dT 




pO 

SS 



r T . r T 



f R ) 



(28) 



The result (J23> is plotted in Fig. as a bold line. This 
indicates that our requirement for \ ^> Fl+Ir results in 
an effectively immediate loss of a portion [fri/(fL-|-fR)] 
of the qubit coherence upon commencing measurement, 
followed by exponential decay at the rate f L . 

It is interesting to note that in the unlikely event of 
the SET island state being known to be occupied at the 
start of the measurement (so that P° — 0), the qubit 
is observed to completely decohere instantaneously. This 
fact alone shows that the quality of charge qubit measure- 
ments by a SET cannot be captured by a simple quantity 
such as rdT m , although this can provide some insight into 
the measurement quality. Note also that g hs( T ) depends 
minimally on the SET response. 



B. Conditional master equation 

We obtain the conditional (or selective, or stochastic) 
evolution of G(t), denoted by a subscript c, in a similar 
way as in Ref. l27l However, we do not eliminate the SET 
island state as in Ref. l27l which was achieved by setting 
r R ,r R > rL,r L . Our result is the following Ito^ con- 
ditional master equation (quantum trajectory/filtering 
equation), where the SET-qubit state is conditioned on 
the tunneling events both on and off the SET island: 



dG c 



y mJ—$ 

^ I Tr [Si 

-dtn hr + iH 



- 1 G 



(29) 



where the non-linear supcropcrator TL is defined by 
TL [x] G ee xG + Gtf -Tr[(x + tf)G] G. The operator 
T and superoperators <Sl an d 5 R are given by 



r; 



h + T L W ® (i - h) 



+ r R S f S <g> h + r R Pb ® (i - h) 



5 L EE T'J 



(£> n 



(1-n) 



T' R J b®n +T K J 6®(i-n) 



(30a) 
(30b) 

(30c) 



In Eq. 1)29(1 we have introduced the classical random 
point processes dNi,(i) and dN R ,(t). They represent the 
number (cither or 1) of electron tunneling events from 
source to island (L) and island to drain (R), respectively, 
in an infinitesimal time interval [t,t + dt). In the sequen- 
tial tunneling mode, at most one of these may be nonzero 
for each particular infinitesimal interval of time. That is, 
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an electron cannot tunnel through the entire SET 'in one 
go', so to speak. The conditional master equation for the 
observable state lf%5)) is identical to (25jl. but with G in 
place of G and H\ in place of H . 

The two terms summed in the first line of Eq. i|29|) 
represent electron tunneling events onto (L) and off (R) 
the SET island. Each event is represented by an incoher- 
ent mixture of the two possible paths by which this can 
occur (corresponding to the two localized electron states 
within the qubit). The final line in Eq. (|29|) represents 
the evolution of G c (t) that is due to the Hamiltonian 
and the null measurement of no SET tunneling events 
in the time interval [t,t + dt). It is important to note 
that these microscopic processes may not necessarily be 
directly observed by a realistic observer 1*22*24 Neverthe- 
less, it is sensible to unravel the master equation lTH)|) 
into a quantum jumpiS stochastic master equation be- 
cause the underlying physical process consists of single 
electron tunneling events. 

The ensemble-average evolution of G(t) [Eq. Ijl9|) ] can 
be recovered by averaging over all possible trajectories 
given by Eq. I|29[) . In the Ito formalism^ this simply in- 
volves replacing the stochastic increments dN^ and dN R 
in Eq. I|29fl with their expectation values, which can self- 
consistently be given in terms of their respective quantum 
averages: 



find from Eq. ljTB|l that 



E[dN hiR (t)]=dtTr[S h , K G c (t)]. 



(31) 



Having presented the quantum trajectory equation 
(|29[1 . we are now in a position to consider the SET mea- 
surement sensitivity using r m and V(t) defined in Sec. 
IIVBI Wc proceed by noting that the qubit state matrix 
p can be represented in terms of the Pauli spin matrices 
as 



j = ^ {P 1 + x a x + y (J y + z a z ) 



(32) 



The conditional master equation [for the observable state 
(I25[) ] can be expressed in terms of the conditional Bloch 
sphere variables x^ , , and well as the conditional 

SET island occupation probabilities P C N . The condition- 
ing here is upon both the measurement results (subscript 
c) and the SET island state (superscript or 1). To find 
V(t) (and hence approximate r m ), we are interested in 
the expectation value of (z c ) 2 = (z° + z\) 2 . Averag- 
ing over the jumps in the conditional master equation 
for z c and using the fact that stochastic variables satisfy 
d(z 2 ) = 2zdz + dzdz, we obtain the following equation 
for the time rate of change of this expectation value: 



dE 



dt 



E 



(2r 



L z c ° + p c °Ar L ) 



4P c °r L + 2Ar L z c ° 

(2f K zl + P^AFr) 
4P c if R + 2Ar R z c i 



21 



(33) 



(AT L ) 2 , ol (AT R y 



4r L 
f L f R 



r 



( AT L 
V2f L 



4r T 



ATf 
2f R 



(34) 



A more precise value for the measurement time is ob- 
tained from Eq. i|13|) by numerical solution of the condi- 
tional master equation for the observable state. Figure^] 
shows plots of V(t) = 1 — E [z 2 (£)] for a strongly coupled 
SET for the same parameter values as the corresponding 
plots in Fig. |3 Notice that the analytical approximation 
(the solid line) agrees very well with all of the numeri- 
cal results (dashed lines). A comparison of Figs. [3] and 
0] shows that the strongly responding SET [in plots (d)] 
has comparable decoherence and measurement times be- 
cause the functions V(t) and g bs(T~) are comparable in 
magnitude at all times. 





Knowing this and using Eq. I|21|l to give zf (0) = 0, we 



t(T[') 



FIG. 4: The variance in a z (the population difference operator 
in the measurement basis), V(t) = 1 — E \z% (t)] . The dashed 
line is an average over 2000 trajectories obtained using the 
conditional master equation. The solid line is the analytical 
approximation assuming exponential decay of V(t). Other 
details are as in Fig. |3| 



C. Quality Of The Measurement 

Previously^iiSiiiiiii^ the quality of the measurement 
of a charge qubit by a (weakly responding) SET has been 
quantified by comparing qubit decoherence and measure- 
ment times, with the quantum limit for a weakly respond- 
ing device given by Eq. QJ. 

More generally, the purity of the conditioned qubit 
state (conditioned by measurement results) provides the 
means to determine the 'ideality' of the detector (how 
close to the quantum limit the detector is operating). 
This is because an "ideal" detector will preserve a pure 
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qubit state throughout the measurement. Thus, we 
maintain that the purity of the conditioned qubit state 
is the ultimate quantifier of measurement quality. The 
conditioned state of the qubit can only be obtained from 
the conditional master equation. 

In this subsection we analyze the quality of the mea- 
surement in two ways. We first use the traditional tech- 
nique of comparing the qubit decoherence (back-action) 
and measurement (sensitivity) times using our definitions 
for these quantities from Sec. I1VI We then analyze the 
average purity of the conditioned qubit state during the 
measurement. 



1. Comparison: Measurement time and Decoherence Time 

Comparing qubit measurement and decoherence time- 
scales gives an intuitive idea of the quality of the measure- 
ment. If these time-scales are of the same order, then we 
can say that the detector is operating close to the quan- 
tum limit. Using Eqs. Ij28(l and Ij34|) we can obtain an 
analytical approximation for TdT m . The result is 



r r iT n 



(i + r L /f R ) 



(Ar L /2f L ) +(Ar R /2f R )' 



(35) 



This result is a valid approximation when V(t) is a de- 
caying exponential, which is valid for a wide range of 
parameter values as shown in Fig. 0] The minimum for 
Eq. (|35|) is 1/2, which occurs only for l\ <C Tr (high 
SET asymmetry) and AFl, R ~ 2rL. R .The second condi- 
tion requires both T' L <C IY and r R <C T r , which could 
in theory be realized by careful arrangement of the SET 
parameters (specifically bias and gate voltages, and tun- 
nel junction conductances). 



2. Average Purity of the Conditioned Qubit State 

Averaging the purity of the conditioned qubit state for 
the four ensembles in Fig. 0] gives the results in Fig. [5] 
Here we choose the qubit to start in a pure state [a super- 
position with z(0) = 0], rather than a mixture. At the 
commencement of the measurement, the purity drops in- 
stantaneously [as did 5obs( T )]- For weaker response [Fig. 
03 (a) and (b)], the purity continues to drop for some 
time because the SET takes much longer to determine 
the qubit state than it does to collapse it. This becomes 
apparent by comparing plots (a) and (b) of Figs. |3|and 
|3J For stronger response the purity approaches 1 rapidly 
[see Fig. (c) and (d)] because the measurement time is 
significantly decreased, as can be seen in Fig. 0] (c) and 
(d). In all four cases shown in Fig. the continuous 
measurement eventually purifies the conditioned quan- 
tum statei 35 ' 36 ! 37 

The competing effects of sensitivity and back-action in 
the qubit measurement can be captured by the purity 
of the conditioned qubit state in the following way. If a 
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FIG. 5: Average purity E \x\ + y\ + zl] over 2000 trajecto- 
ries. The four plots correspond to the same ensembles calcu- 
lated in Fig. 0] 



pure qubit state remains pure during the measurement, 
the detector is operating at the quantum limit and can be 
called an ideali detector. In order for the qubit to remain 
in a pure state once the measurement starts, the SET 
island must be unoccupied when in steady-state. This 
occurs when electrons tunnel off the island almost imme- 
diately, i.e. when r^T^ < r R ,r R . This "asymmetry" 
adiabatically eliminates the influence of the SET island 
state on the qubit, as in the model of Ref. lU With this 
condition satisfied, we obtained another ensemble of 2000 
quantum trajectories using the conditional master equa- 
tion and plotted the average purity of the conditioned 
qubit state in Fig. El Note the different y-axis limits. 
Again the SET response was controlled by changing 
[again it was decreased incrementally from (a) to (d)]. 
The SET response was, from (a) to (d), AI/I = 0.23, 
0.58, 1.32, and 1.70. Other parameters for Fig. were 

n = o.oir L , r R = ioor L , r R = ioor R , x = 5or R . 

Similar results were obtained for T R = 10r R . For the 
strong response plot of Fig. (d), this set of parameters 
corresponds to the regime given after Eq. (|35|l . That is, 
the regime where the quantity TdT m approaches 1/2. In- 
deed, the values of FdT m for Fig. [B] approach this limit of 
1/2, as given in the figure caption. 

Figure [5] shows that the average purity of the condi- 
tioned qubit state remains very close to 1 throughout the 
measurement for the asymmetric SET. Thus we can say 
that the strongly coupled asymmetric SET operates very 
close to the quantum limit. This conclusion could not 
be drawn from the values of FdT m alone. Interestingly, 
the SET response [increasing from plot (a) to (d), as be- 
fore] has very little effect on a SET with high asymmetry. 
This is similar to the case of a QPC, which is an ideal 
detector independent of the strength of its response. To 
summarize, the conditions required for a strongly cou- 
pled SET to operate close to the quantum limit are: (a) 
< < T L , and (b) f L < f R . 
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FIG. 6: Average purity of the conditioned qubit state, 
E [a;;? + Vc + %c] , when monitored by an asymmetric SET 
(Fl <C Tr). The average is over 2000 trajectories. Param- 
eter values are given in the text. From (a) to (d), we have 
r d r m = 1.03, 0.95, 0.72, and 0.59, respectively. 



VI. CONCLUSION 

In this paper we have examined sensitivity and back- 
action in charge qubit measurements by the single- 
electron transistor (SET) operating in the sequential tun- 
neling mode. We have shown that the unconditional mas- 
ter equations derived in previous work2ii2ii2ii£ii& for this 
measurement (which are equivalent to each other) are 
not of the Lindblad formii that is necessary and suffi- 
cient for guaranteeing valid physical states. Thus, these 
master equations cannot be used to model real physical 
systems. In particular, this statement also applies to the 
equation for the conditioned state (one that an experi- 
menter would calculate knowing the output of the detec- 
tor) in Ref. Il^2i The unconditional master equation was 
derived independent of the strength of SET-qubit cou- 
pling, but analysis of the SET measurement quality (or 
efficiency) assumec&iSiiiii^ weak SET response. Weak 
response is when the average component of the SET out- 
put is much larger than the component that depends on 
the qubit state. 

We have presented an alternative master equation that 
is of the Lindblad form, as well as the associated quan- 
tum trajectory equationi£*i2*22i for the conditioned state. 
To derive these equations required assuming strong SET- 
qubit coupling. A strongly coupled detector can exhibit 
strong response or weak response, whereas weak coupling 
implies weak response. 

We have shown that, in general, quantifying the qual- 
ity of charge qubit measurements by a SET is not as 
simple as merely comparing the qubit decoherence (also 
called coherence) time and measurement time (as done 
previously for the weakly responding SET in Refs. 0, 
EUll and This is because, for strong SET-qubit 
coupling, the SET island charge state has a large effect 
on the qubit coherence. Despite this, comparing qubit 



decoherence and measurement times does provide some 
insight into the measurement quality. 

Previous definitions for the qubit decoherence 
ti ma &i2iiiLA2ii£ either produced an overestimate of the 
measurement quality (as in Refs. 0,0,0, andll5T). or 
are only valid for a weakly responding detector (as in 
Ref. We presented a new definition for the deco- 

herence time that is valid for strong SET-qubit coupling 
and arbitrary response, by adopting a coherence function 
technique^ from the field of quantum optics. 

Previous definitions for the measurement 
timeiiiliidtli are only valid for weakly respond- 
ing detectors. Thus, we also presented a new definition 
of the measurement time which is analogous to the 
new decoherence time definition. Our definition is in 
terms of the conditioned qubit state, and involves an 
average over many stochastic measurement records. 
An analytical approximation for the measurement time 
compared well to the numerical solution. A similar 
approach to quantifying the measurement sensitivity 
was first considered in Ref. I27I 

An important result of this paper was to show that the 
strongly coupled SET can operate at, or very near, the 
quantum limit during charge qubit measurements. This 
was shown in two ways. The first involved comparing 
the decoherence and measurement times found using our 
Lindblad-form master equations. In the strong response 
limit, these times were found to be of the same order as 
each other. For the highly asymmetric SET where tun- 
neling into the drain (collector) occurs at a much higher 
rate than tunneling onto the SET from the source (emit- 
ter), these times were of the same order as each other, 
independent of the SET response. This is also the case for 
the QPC, which can be operated as an ideal detector^ 
The second, and our preferred, way to show that the 
strongly coupled SET can operate close to the quantum 
limit was by considering the average qubit purity for an 
ensemble of stochastic measurement records. Again the 
asymmetric SET displayed near ideal detector behavior 
by maintaining the purity of the qubit state above 98%. 

As far as we are aware, our Lindblad-form quantum 
trajectory equation is the first of its type presented for 
charge qubit measurement by a SET. Moreover, having 
shown that the SET can approach the quantum limit for 
charge qubit measurements, it is now meaningful to study 
how the measurement is affected by extra classical noises 
and filtering due to a realistic measurement circuit 
This research is continuing. 
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APPENDIX A: SENSITIVITY AND 
BACK-ACTION IN A QPC MEASUREMENT 

The qubit decoherence rate for measurement by a low- 
transparency QPC (single tunnel junction device) can be 
obtained either by inspection of the master cquation^2ii£ 
or using the technique introduced in Sec. II V Al The same 
result is obtained using both methods: 



(Al) 



where V (r') is the tunneling rate through the QPC when 
the target dot is unoccupied (occupied). Using the new 
technique introduced in Sec. UVBl the analytical approx- 
imation of the qubit measurement time (which agrees 
with Ref.OU is 



|AT|- 

4r 



Vr- 



r + r 



Combining these results gives 

r + r 



r r iT 1T 



Vr- 



2 ■ 



(A2) 



(A3) 



Note that 1 > TaT m > 1/2. The upper bound is ap- 
proached for a "strongly responding" QPC (r ^> V). 



The lower bound is approached for a "weakly respond- 
ing" QPC (r w r'). i.e. when the measurement and 
decoherence rates go to zero. 



APPENDIX B: THE SET-QUBIT MASTER 
EQUATION OF GOAN 

In Ref. ITri Goan presents a master equation for the 
combined state of the qubit plus island and drain of the 
SET. In that paper, the master equation is shown to be 
equivalent to that in Ref. 0. Using our notation and 
tracing over the drain degrees of freedom gives 



G = 



H,G 



T' L V 



G + T'V 



6 <g> n 



S f <g> (1 - n) 



G + TrD 



6® (1 - h) 



+2r L D [h] PGb + 2T R V [h] bGP, 



G 
G 

(Bl) 



where V 1 corresponds to e 1 ^ in Ref. ITfil and we have used 
[h, [h, x}) = —22? [n] x. This master e quat ion is equiva- 
lent to those in Rcfs. H E3, Q and M Comparison 
of this equation with our unconditional master equation 
(|19fl shows that the terms in the final line of Eq. I)B1|) 
are not present in our model. These are the terms that 
lead to the non-Lindbladian dynamics in Refs. I3L IT2I and 
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